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We probe for statistical and Coulomb induced spin textures among the low-lying states of 
repulsively-interacting particles confined to potentials that are both rotationally and time-reversal 
invariant. In particular, we focus on two-dimensional quantum dots and employ configuration- 
interaction techniques to directly compute the correlated many-body eigenstates of the system. We 
produce spatial maps of the single-particle charge and spin density and verify the annular structure 
of the charge density and the rotational invariance of the spin field. We further compute two-point 
spin correlations to determine the correlated structure of a single component of the spin vector field. 
In addition, we compute three-point spin correlation functions to uncover chiral structures. We 
present evidence for both chiral and quasi-topological spin textures within energetically degenerate 
subspaces in the three- and four-particle system. 
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I. INTRODUCTION 

The investigation of correlations among electrons con- 
fined to quantum dots (QDs) is an active area of research 
in conden sed matter physics due to their experimental 
tunabilityP^ their theoretical emcacyP and their appli- 
cation in, for example, quantum information scienceP"21 
For two-dimensional QDs with circular confinement, the 
combination of confinement and long-range Coulomb re- 
pulsion results in charge densities peaked in annular re- 
gions about the dot center. Calculations of two-point cor- 
relation functions^ have further revealed that these con- 
fined electrons exhibit textures akin to Wigner molecules. 
Numerical work has shown these behaviors to depend on 
the size and shape of the quantum dot, as well as the 
strength of the applied magnetic fieldP^ 

Since the spin of confi ned el ectrons provides a viable 
implementation of qubitsp^ an understanding of the 
configurations and correlations formed among confined 
spins is crucial for the implementation of spin-qubits. 
Furthermore, the degree of control in fabrication and ma- 
nipulation of QDs makes them ideal environments for the 
study of fundamental behaviors of both spin and charge. 
In this paper, we investigate the correlations that exist 
between the spins of electrons trapped in circular QDs. 
We are specifically interested in the formation of topo- 
logical spin textures that may arise due to interaction or 
statistical effects among the confined charges. 

A significant challenge in the development of spin- 
based quantum dot quantum computing is the suppres- 
sion of decoherence of the spin-states for time-scales 
much longe r than the time required to controllably flip 
a spinp3LL21 The potential to encode information using 
topological degrees of freedom is appealing since the 
enhanced stability can mitigate the burden of error- 
correction. Various schemes which exploit a system's 
topological structure have been proposedP^S^ The ad- 
vantage of these is that they are physically fault-tolerant; 



they are immune to local perturbations that degrade the 
coherent evolution of the state, a necessary ingredient in 
quantum computationPHl Possible systems include two- 
dimensional spin modeli^ (for example, atoms in an op- 
tical lattice^) and fractional quantum Hall systems.^ 
These proposals rely on the existence of non-Abelian 
anyons in the excitation spectrum of the models for the 
information processing. 

However, even in the absence of anyonic excitations, 
textures with topological structure are expected to be 
long-lived in-and-of-themselves due to their global cor- 
relations. Since correlations decrease decoherencej22H2l] 
these topological structures could form an important pro- 
cessing element in more conventional quantum comput- 
ing schemes. Even in finite-sized systems, where true 
topological stability likely does not occur, the relevant 
relaxation and decoherence times can be significantly en- 
hanced. In Ref. [22j for example, even moderate charge 
correlations were sufficient to more than double the de- 
coherence time. 

Numerical work has predicted the formation of spin 
textures in QDs immersed in a magnetic fieldPS Ex- 
perimental evidence suggests the existence of fermionic 
spin textures in a two-dimensional electron gas (2DEG) 
confined in semiconductor heterostructures, 26 in vertical 
QDsJ^l and in few-electron lateral QDs! 2 ^ 

A topological spin qubit would be advantageous as it 
could be more robust against local environmental per- 
turbations. Nuclear magnetic resonance measurements 
of GaAs/AlGaAs quantum wells have shown evidence for 
the localization of topological sk yrmio n spin textures as 
the temperature approaches K) 30 * 31 ! Recently, topolog- 
ical spin textures have been experimentally observed in 
topological insulators by means of spin-resolved angle- 
resolved photoemission spectroscopy.^ The emergence 
of a skyrmion lattice has also been detected in the chi- 
ral magnet MnSi using neutron scattering.^ Topologi- 
cal textures that are predicted to appear in QDs include 
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vortice^HEH anc j merons l31H2I] Vortices occur in the pres- 
ence of a strong external magnetic field, when electron 
current circulates in a plane around localized regions of 
low electron density. Merons are topological spin textures 
characterized by a central "up" or "down" spin which 
smoothly transitions into an in-plane 27r-winding along 
its boundary.^ As developed in the theory, the realiza- 
tion of both types of quasiparticles requires the presence 
of an external magnetic field. 

In this work, we present evidence for the existence of 
spin textures in circular QDs for both three-electron and 
four-electron systems in the absence of an external mag- 
netic field. The electronic wave functions are calculated 
by configuration interaction techniques. Two-point and 
three-point spin correlations are calculated in order to 
uncover both correlation and chirality in the spin tex- 
tures which are concealed in superpositions of different 
configurations. 

In Sec. [Hj we introduce our model and the foundation 
upon which our calculations are based. Section III de- 



scribes specifically the spin correlation calculations used 
to examine the spin textures in the QD. We then go on 



to describe our results for systems of three (Sec. IV I and 



four (Sec. [VJ) interacting particles. We conclude with a 
summary of our findings, their implications, and sugges- 
tions for further investigations in Sec. |VI| 



II. QUANTUM DOT SYSTEM 

Our system consists of N interacting quasiparticles of 
charge e, bound to a two-dimensional (2D) plane and lat- 
erally confined by a parabolic potential. The 2D Hamil- 
tonian used to describe this "standard model" is 
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where e is the dielectric constant of the medium and h is 
the single-particle Hamiltonian describing harmonic con- 
finement; 
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where m* is the effective mass, r = (x, y) the position 
operator, and luq the parabolic confinement frequency. 
Throughout this paper, we take the magnetic field to be 
zero, and therefore set the vector potential A (f ) = 0. 

Two harmonic-oscillator quantum numbers, n, m = 
0,1,2,..., characterize the eigenstates of the single- 
particle Hamiltonian, 3?' Eq. (lb I. These eigenstates are 



the "atomic orbitals" of the QD, and are given by 
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where, a? and are the usual Bose creation operators, 
and 1 00) is the single-particle ground state. These or- 



bitals have energy e nm given by 

£„ m = 7jfi + (n+A)+7LO_(m+ A), (3) 



where ft± — (\/4wq + u> 2 ±w c )/2, and uj c = eB/(m*c) is 
the cyclotron frequency. This energy reduces to hu)o(n + 
m + 1) in the absence of a magnetic field. 

The single-particle Hamiltonian, the z-component of 
the orbital angular momentum, L z , and a component of 
the spin operator — which we take to be the ^-component 
S z — form a set of commuting observables which we take 
to classify our states: L z \nms) — h{n — m)\nms), 
S z \nms) = hs\nms), h\nms) = e nm \nms). 

We are interested in spatial textures formed by the spin 
field, and so we require the position-space representation 
of the orbitals.^ These are given bj@3l 



<f>nm(r,e) = (-1) 
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where r and are the polar coordinates in two dimen- 
sions, Iq — (2m*uj) is the effective length with uj = 
\J 4ujq + w 2 /2, n! — n+m, mf — n—m, n r — (p! — \m'\)/2, 

and Ln (x) is the generalized Laguerre polynomial.^ 

The eigenstates of the interacting system are deter- 
mined by exact diagonalization of Eq. (la). This pro- 



cedure begins by determining many-particle basis states 
(Slater determinants) , eigenstates of Eq. ( lb ) , that are 



composed of antisymmetrized products of the single- 
particle states in Eq. pi). We use 288 single-particle 
states, and as many as 4900 many-particle basis states 
in the diagonalization routine. Without loss of fidelity, 
and for computational efficiency, the number of many- 
particle basis states is reduced when determining the 
two-point and three-point spin correlation calculations 
over the range of the entire QD. Block-diagonalization 
is performed for a given set of parameters. These in- 
clude system parameters (B, uj q , to*, e) and the con- 
served quantities N, L z , S z , S 2 . The Coulomb matrix 
elements are evaluated using the convenient closed-form 
expression derived in Ref. 1461 

With the eigenstates determined, we calculate one-, 
two-, and three-point position-dependent spin correlation 
functions over energetically-degenerate manifolds. The 
structure of the particular operators used in these calcu- 
lations is discussed in the next section. 



III. PRODUCT SPIN OPERATORS 

For our investigation, we require the products of up 
to three one-body spin operators. Here, we introduce 
the product-operators used in the calculations shown in 
the proceeding sections. The details associated with the 
derivation of each product-operator are discussed in the 
Appendix. 
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A. One-Body Spin Operators 



B. Two-Body Spin Operators 



Except where indicated, all averages are taken over 
energetically-degenerate manifolds. That is to say, 
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where the density operator p is defined as 

n 
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and where the states \Ei) are all the states in a given 
degenerate manifold of Eq. |l]): H\Ei) = const. \Ei) for 
all? = 1, 2, . . . n. 

Our analysis begins with the evaluation of both spin 
density and number density at position r in the QD sys- 
tem. We are specifically interested in isolating the spin- 
up and spin-down densities along the coordinate axes. 
Due to the conservation of spin in the system, these one- 
body spin operators can only distinguish the spin-up den- 
sity from the spin-down density along a single axis. We 
therefore define a set of spin operators S± z (r) that sep- 
arately determines the spin-up and spin-down densities 
along z at position r. In canonical form (see Appendix 
Eq. ( Al )), this is given by 



t (r) = TU 



(6) 



where U+. 



scj)*{r)(i) j (r)/2, with s = ±1 (ft = 1). Note 



the composite indexes i and j each represent a set of 
orbital quantum numbers n and m. 

The number density and spin density operators at po- 
sition r are then given by 

n(r) = (r) = 2 (S +Z (r) + S- Z (r)) , (7) 

S 

Sz{r) = ^$>$W» = (S +Z (r) ~ S- Z (r)) , (8) 

s 

where tpl(r) and ip s (r) are field operators that respec- 
tively create and annihilate a fermion at position r, with 



spin s. In terms of the eigenstates \nms) of Eq. ( lb I, the 
field operators are given by 
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with 4> nm (r) = (r\nm) given in Eq. Q. 

The effects of Coulomb interaction between the par- 
ticles are apparent when the expectation values of the 
above operators are compared between the interacting 
and noninteracting systems. 



We next investigate the two-point correlations as pro- 
jected onto the z-axis. (Any other choice yields iden- 
tical results.) Unless otherwise indicated, a spin-up 
(spin-down) particle refers to a particle with spin po- 
larized along the positive (negative) axis of quantization 
(in this work, the z-axis). Specifically, we investigate 
(S +z (r )S +z (r 1 )) and (S +x (r )S- z (ri)). In canonical 
form, the operators are given by 



S +z (r )S az ( ri ) = 
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with 



Qij'u = ^ (0i(ro)^(r 1 )^ fc (r o )^(n)M-'«r 

and a = ±1. As in Eq. ([6]), the indexes i through I in 
Eq. ( 10 1 are again composite indexes of pairs of orbital 



quantum numbers. The two-point spin correlations mea- 
sure the probability of finding a particle with spin pro- 
jection az at position r\ given the existence of a particle 
with spin projection +z at position Tq. 

Finally, the correlation between a spin-up particle at 
r and the net spin density at ri is 

S +z (r )S nct z (ri) 

= S +z (r )S +z (r{) - S +z (r )S- z (n). (11) 

These two-point spin correlations are useful for deter- 
mining parallel or antiparallel spin properties such as 
magnetic ordering.^ They are insufficient, however, for 
determining chiral textures where correlations are mea- 
sured with respect to orthogonal axes. For that we turn 
to the three-point spin correlation functions next. 



C. Three-Body Spin Operators 



We compute three unique three-point correlations: 



and 



(S +x (r )S +y {r 1 )S nct x (r 2 )) , 
(S+xi^S+yir^Snct v (r 2 )) , 



(S +x (r )S +y (n)S net z (r 2 )) 



(12a) 
(12b) 



(12c) 



These three-point spin correlations measure the proba- 
bility of finding a particle with spin projected along the 
x, y or z axis, respectively, given there is a particle at po- 
sition ro that is spin-up along the x-axis and a particle 
at position r\ that is spin-up along the y-axis. Whereas 
the two-point functions can determine whether the spin- 
projection of a second particle is parallel or antiparallel 
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to the spin-projection of the first particle, it cannot de- 
termine the orientation of the spin of the second particle 
in a plane other than that of the spin of the first par- 
ticle. The three-point functions in Eq. ( 12 ) can indeed 



uncover such chiral structure. Explicitly, the three-body 
spin operators can be expressed as 

net x (ra) 

— — ?T s i s 2 s 3-t ~t at A A {-[la} 
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x [is 3 (/)*(ro)(t)*(ri)(t)l(r 2 )(t) m (ro)(l)n(ri)(l)i(r2) 
- ^ l *(T , o)^fe(r"i)^(r"2)0m(ro)0;(ri)0„(r 2 ) 

-^(r^ir^ir^MMr^M^)] , (13b) 

and with similar expressions for the remaining two op- 
erators in Eq. (12). In Eq. (13 1, each of the six indexes 



i through n is once again a composite index over pairs 
of orbital quantum numbers. Furthermore, in Eq. (13b I, 
Si = — S{, with Si = ±1 denoting the usual spin projec- 
tions. 

The choice for the first two spin-projections is not 
unique; due to the absence of a preferred spin orienta- 
tion in the system, each expression is equivalent to any 
cyclic permutation of the spin components. We focus be- 
low on the cases where two of the three spins operators 
lie in the (x-y) plane of the dot. 



IV. THREE-PARTICLE SYSTEM 

In this section, we investigate spin correlations that ex- 
ist in the two lowest-lying states of a system with three 
charged particles. Our system is modeled with GaAs 
parameters (e = 12.4 and m* = 0.067m e ), and our con- 
finement potential is wo = 1.0 meV, yielding an effective 
length at B = of l = 23.8 nm. 



A. Ground State Manifold 

At zero magnetic field, the ground state of the three- 
particle system is four-fold degenerate, with quantum 
numbers L z = ±1, S = 1/2, and S z = ±1/2. We 
compute (S +Z (r)) and (S- Z (r)) within this degenerate 
subspace. [See Eq. From these, we determine both 
the net density, (n(r)), Eq. and the net spin (S z (r)}, 
Eq. We then go on to calculate the two-point spin 
functions to demonstrate correlations between parallel 
and antiparallel spin components, followed by the three- 
point functions to uncover chiral correlations. 
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FIG. 1: (Color online.) Radial dependence of the single- 
particle spin and number density in the three-particle ground 
state manifold. Inset shows densities in absence of Coulomb 
interaction. (Inset axes have the same dimensions as the main 
plot). 



1. Single-particle densities 

To illustrate the effects of long-range Coulomb repul- 
sions, we consider the spin density with and without in- 
teractions. In the non-interacting case, each eigenstate 
is a single antisymmetrized orbital configuration. For 
the three-particle system, there are two particles on the 
\nm) = 1 00) orbital, and one on either the 1 10} or 1 01) 
orbital. This yields four degenerate states with quantum 
numbers (S,S Z ,L Z ) = (1/2, ±1/2, ±1). For the inter- 
acting case, these symmetries are not explicitly broken; 
the degeneracy and the quantum numbers remain the 
same, but the states themselves are now correlated, in- 
volving many other orbital configurations consistent with 
the symmetry. 

Figure [T] shows single-particle densities for the to- 
tal ground-state manifold as a function of radial dis- 
tance from the center for both the interacting and non- 
interacting cases. There is azimuthal symmetry for these 
configurations due to the underlying circular symmetry 
of the dot itself, manifest in the Hamiltonian, Eq. (JlJ. 

The non-interacting case is characterized as Gaussian- 
like with a peak at the origin. When Coulomb interac- 
tions are considered, the repulsion smears out the density 
over different orbital configurations.^ The competition 
between repulsion and confinement results in an annular 
density about the origin. These interaction effects are 
strong; the ground-state energy of the interacting system 
is 10.30 meV for these experimentally-relevant system 
parameters — more than twice the ground-state energy of 
the non-interacting case (4.0 meV). 

The effects of Coulomb repulsion are also reflected in 
the single-particle spin densities as well. Note, how- 
ever, that in both cases we have (S +Z (r)) = (S- Z (r)): 
(Snot z{ r )) = everywhere in the dot, therefore the 
one-point calculations are insufficient for showing the 
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(S +z (r )S +z (ri)) 
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FIG. 2: (Color online.) Two-point spin correlation functions 
for the three-particle ground state manifold. Distributions 
are shown for the spin-up (a) and spin-down (b) densities at 
ri given a spin-up particle at ro = (39 nm, 0), denoted by a 
cross. 
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FIG. 3: (Color online.) Trace along r max = 39 nm in the 
three-particle ground state manifold, revealing the spin dis- 
tribution with respect to a spin-up particle at ro = (r max , 0). 



Coulomb effects on spin. This is a consequence of the 
SU(2) symmetry present at zero magnetic field. 



2. Two-point spin correlations 



For the two-point spin correlations, Eq. (10 1, we con- 



sider the case where ro is fixed at the location of maxi- 
mum single-particle density, r max = 39 nm, as obtained 
in the previous section. We further define the angular 
location of r to be 8 = 0. We do not the discuss the 
non-interacting limit for these calculations. 

The two-point correlations (S +z (r )S +z (ri)) and 
(S+ x (ro)S- x (rx)) for the interacting ground-state man- 
ifold are shown in Fig. [2] as a function of r\ for ro = 
(Vmaxj^o)- Two peaks are evident along the ring of ra- 
dius r max . Note that our averages, Eq. (5), are obtained 
by tracing over all the degenerate states in the ground- 
state manifold. An incipient Wigner crystallization is 
apparent with the spins forming a classical-like lattice at 
the vertexes of a triangle? 10 * 11 ! This structure is not seen 
in the non-interacting case, implying that the Coulomb 
repulsion between the particles is responsible for this spin 
texture. 

To more clearly probe the angular inhomogeneity, we 
plot in Fig. [3] results along the ring r max . In particu- 
lar, wc show the net spin S z as well as the individual 
components S± z , given a spin-up particle at ro. Note 
as ri approaches ro the remaining spin- up density goes 
to zero, indicative of a Pauli vortexPS at that position. 
As well, the spin density at the two peaks is not fully 
polarized, indicating a degree of canting away from the 
z-axis: The net spin tilts towards the x-y plane. The lack 
of equal magnitudes of spin-up and spin-down probabil- 
ities at every point along r max in Fig. [3] indicates that 
the spin density never lies completely in the x-y plane. 
Since the spin-density never crosses through the plane, 
it cannot have winding order. Windings are important 
in these spin systems as they represent clear examples of 
topologically stable structures. 

The two-point correlations are insufficient to determine 
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FIG. 4: (Color online.) Canting angle a with respect to the 
axis of quantization (the z axis). The spin is not fully polar- 
ized along either the +z or —z axis but has components of 
both projections, thus smearing out along a cone of probabil- 
ity about the z-axis. 



the probable orientation of the spin in the x-y plane, so 
the results of these calculations can be interpreted as the 
smearing out of the net spin across the surface of a cone 
centred at a local z-axis at every point measured in the 
calculation. The opening angle of the cone is twice that 
of the local canting angle, the degree of tilting from the 
local z-axis, of the spin. This cone, along with the canting 
angle, is illustrated in Fig. |4j 

The canting angle can be determined in the following 
manner: If we consider the spin density field shown in 
Fig. [3] as itself a spin-half field, we may generally write 
its local orientation as 

Mw.fli)) = i (c + \+z)+e^c.\-z)) , (14) 



where c± are real and may be defined as 



c+ = (S +z (r )S +z (r n 



l)) = iVcos 



(15a) 



c_ = (S +2 (r o )£U(r max ,0i)> = N sin-, (15b) 



and where N 



c~_ is a local normalization. The 
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FIG. 5: (Color online.) Canting angle a as a function of 6\ 
along the ring of maximum single-particle density r max = 39 
nm in the three-particle ground state manifold. The net spin 
density from the two-point calculation along the same ring 
r m „ is shown for reference. 
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FIG. 6: (Color online.) Three-point spin correlations, 
Eq. (12) , along the ring r max = 39 nm for the three-particle 
ground state manifold, (a) The net S y is plotted given a spin- 
up particle along the a>axis at 6q — 0° and a spin-up particle 
along the j/-axis at 9\, (b) The three-point correlations for 
net spins along the as, y and z axes for Oo — 0° and 9\ = 126° 
on the ring r max . 



symmetry of the two-point functions prevents discrimi- 
nation of different values of the azimuthal angle /?, but 
it can determine the canting angle a. 

For each of the two peaks in Fig. [3j the canting angle 
with respect to the positive z-axis is determined to be 
a = 131°. That is, relative to a spin- up particle at r , 
the spin-density peaks describing the other two particles 
both occupy the surface of a cone with canting angle 49° 
from the negative z-axis. Figure [5] shows the resulting 
canting angle from the positive z-axis of the local Bloch 
vector for each value of 0\ along the ring r max . The 
canting angle becomes greatest (a approaches 180°) as 
6*i approaches 0q = 0, and is a minimum at 0-y = 180°. 
(Pauli exclusion dictates that a —> 180° as 0\ — > 0q.) 

The two-point functions in this system with SU(2) 
symmetry are insufficient to distinguish chiral structures. 
Three-point functions are necessary to resolve spin com- 
ponents in the plane perpendicular to the axis defined by 
the two-point functions. We turn to these next. 




FIG. 7: (Color online.) Net spin distribution in the plane of 
the QD for the three-particle ground state manifold, as de- 
termined by three-point spin correlations. The magnitude is 
denoted by the color bar and the direction is indicated by the 
vector field. The ^-component is negligible here. The calcula- 
tion is done with respect to a particle of spin-up projection on 
the a>axis at r$ — (39 nm, 0°) and a particle of spin- up pro- 
jection on the i/-axis at ri = (39 nm, 126°). (These locations 
are denoted by the blue arrows. 



3. Three-point spin correlations 



As described in Sec. IIIC we compute the three dis- 
tinct three-point correlation functions given in Eq. (|12|. 



Other three-point functions can be related to one of these 
three due to the symmetry of the underlying Hamilto- 
nian. 



In Eq. ( 12 ), we fix r$ on the ring of maximum single- 
particle density (r max = 39 nm) at 9 = 0. The cor- 
relation functions are then a function of the four vari- 
ables r\, 0i, r%, and 02- If we further choose to 
probe the system along the ring r max , we then obtain 
the two-dimensional map in the two angles 0% and 0i 
shown in Fig. |6^a). Explicitly, Fig. [6]ja) is a plot of 
0)S +y ( 

^maxi ^l)^nct yf^maxj 62)) as a function 
of 0\ and 02- As the angular location of the second spin 
approaches 0\ — 120° with respect to the first spin, a 
peak emerges in the net spin distribution along the y- 
axis. The maximum correlations occur at (0\,02) = 



(126°, 235°) and (234°, 125°). Note the inversion sym- 
metry about the point (180°, 180°). Very similar results 
are seen for the net spin distribution along the x-axis 
(not shown). 

The three-point correlations for the net spin distribu- 
tions along the x, y and z axes with respect to a spin-up 
particle along the a;-axis at r$ — (r max , 0), and a spin-up 
particle along the y-axis at r\ = (r max , 126°) are dis- 
played in Fig. [6jb) . Note the net spin distribution along 
the z-axis is negligible in comparison to the distributions 
along the x and y axes: Any spin density that remains 
in the system lies only in the plane of the two spins at r 
and n, respectively. 

Figure [7] is a map of the net spin distribution in the 
three-particle ground-state manifold, given a spin-up par- 
ticle along x at r$ — (r max , 0) and a spin-up particle 
along y at T\ = (r max ,I26°). The peak along the ring 
r max occurs at #3 = 243°, and the net spin at this peak 
has an equal spin-down projection along each the x and 



7 



y axes. Note in this distribution the peak of the net spin 
density and the two locations ro and r\ are each equidis- 
tant from each other. 

The three-point correlations suggest that the most 
probable spin configuration in this three-particle state 
has a planar, splayed order, as seen in Fig. [7] Recall, 
however, that the results from the two-point correlations 
imply that the state does not have a clear winding order. 



B. First Excited State 

Like the three-particle ground state, the first excited 
state is also four-fold degenerate. However, in contrast 
to the ground state manifold, L z = for all four de- 
generate states in the first excited state. Degeneracy oc- 
curs through the spin quantum numbers S — 3/2, and 
S z = ±3/2, ±1/2. 

As in the ground state, the single-particle density and 
the spin density distributions are rotationally symmetric. 
When interactions are considered, the annular distribu- 
tions peak at r max = 41 nm. These distributions are 
similar to those seen in Fig. [I] for the ground state and 
are not shown. As seen in the ground state manifold, 
the net S z is once again zero everywhere in the QD for 
both the interacting and non-interacting cases. The non- 
interacting limit is not investigated further. 

The two-point spin correlations in the interacting first 
excited state are similar to those seen for the interacting 
ground state (see Fig. |2j, namely, there are two peaks 
along the ring r max . The two-point distribution also 
shows canting towards the x-y plane with respect to a 
spin-up particle at ro = (r max , 0). Unlike the ground 
state, this canting is constant for all 9 along the ring 
r max except at a very small region about 9 — in accor- 
dance with Pauli exclusion, and has a value of a = 53°. 
Once again a Pauli vortex is present at r . 

We turn now briefly to the three-point spin correlations 
in the first excited state along the ring r max . We fix one 
spin-up particle along the z-axis at 9q = 0°, and the 
other spin- up particle along the y-axis at 9\ = 120°, the 
location of one of the two peaks determined by the two- 
point correlation calculation. Due to the spin polarized 
states, minor differences exist between this manifold and 
the ground state manifold, but otherwise these results are 
very similar and so results for the excited state manifold 
are not shown. 



V. FOUR-PARTICLE SYSTEM 

We examine the lowest two energy eigenstates in a sys- 
tem of four charged particles for the same system param- 
eters used above. We begin the investigation of the four- 
particle system by first determining the radial location of 
maximum single-particle density in each state. With this 
established, we calculate higher-order spin correlations in 
this region. 
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FIG. 8: (Color online.) Two-point spin correlations in the 
four-particle ground state manifold given a spin-up particle at 
ro = (48 nm, 0) (white cross). The distributions are shown 
for the net spin- up density (a) and spin-down density (b). 



A. Ground State 

The ground state of the four-particle system is three- 
fold degenerate with quantum numbers L z = 0, S = 1, 
and S z = ±1; 0. The single-particle density distribution 
in the ground state manifold (not shown) is annular in 
shape with a peak at r max = 48 nm from the center of the 
dot. This distance is greater than in the three-particle 
states primarily due to the additional Coulomb repulsion 
present in the system. 

The two-point spin correlations in the four-particle 
ground state are calculated with respect to a spin-up par- 
ticle at ro = (r max , 0). The calculations are similar to 
those for the three-particle system. Figure [8] shows the 
distribution of the remaining spin-up density and spin- 
down density in the ground state manifold with respect 
to the spin- up particle at ro. As seen in the three-particle 
ground state, there is a Pauli vortex at r . Figure [8]( a) 
shows that the probability of finding another spin-up par- 
ticle is strongest at 9\ — 180°. Conversely, the spin- 
down distribution in Fig. [8jb) shows two peaks, one at 
9 1 = 90° and the other at 6>i = 270°. The saddle point 
at 9\ = 180° is more than half of the magnitude of the 
peaks. Taken together, the two plots in Fig. [8] indicate an 
antiferromagnetic alignment of the spins, with each spin 
equidistant from each other, distributed along the ring of 
maximum single-particle density. This is further revealed 
in the net spin distribution, shown in Fig. [9} Figure [9^b) 
shows the evident antiferromagnetic tendency in the four- 
particle ground-state manifold. However, as shown in 
Fig. 9ja), in this small system, the spins are neither fully 
localized nor fully polarized in the quantum dot. 

In the four-particle ground-state manifold, the net spin 
density from the two-point calculation shows canting 
along the ring r max as a function of 9\ similar to that 
in the three-particle ground-state manifold. The cant- 
ing angle becomes antiparallel (a — 180°) relative to the 
spin-up particle at ro — (r max , 0) as 9\ approaches zero, 
passing through the x-y plane when the net spin den- 
sity from the two-point function is zero along the ring 
r max (near 9 X = 130° and 9 X = 230°), and has a positive 
projection on the z-axis at 9\ = 180°, where the two- 
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FIG. 9: (Color online.) Two-point spin correlations in the 
four-particle ground state manifold, (a) Trace is along the 
ring of maximum single-particle density r max = 48 nm, with 
respect to a spin-up particle located at ro — (r max ,0°). (b) 
Net spin in the QD, given a spin-up particle at ro = (r max , 0°). 




120 180 240 
0i (degrees) 



FIG. 10: (Color online.) Canting angle a as a function of 6i 
along the ring of maximum single-particle density r max = 48 
nm in the four-particle ground state manifold. The net spin 
density from the two-point calculation along the same ring 
r max is shown for reference. 
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FIG. 11: (Color online.) Three-point spin correlations along 
the ring of maximum single-particle density r max = 48 nm for 
the four-particle ground state. Specifically, the net S x (a) and 
net S y (b) are determined with respect to a particle of spin-up 
projection along the x-eotis at ro = (r max , 0) and a particle 
of spin-up projection along the y-axis at n = (r max ,#i). 
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FIG. 12: (Color online.) Three-point spin correlations along 
the ring of maximum single-particle density r max = 48 nm 
in the four-particle ground-state manifold with respect to a 
particle that is spin-up along the x-axis at ro = (r max , 0) and 
a particle that is spin- up along the i/-axis at ri = (r max , 90°) 
(a) and n = (r max , 180°) (b). 



point function shows the net spin density to be predom- 
inantly spin-up. The canting angle along r max is consis- 
tent with the antiferromagnetic ordering suggested by the 
two-point correlation calculations. The result is shown in 
Fig. [101 

To investigate possible chiral textures, we examine 
three-point spin correlations in the four-particle ground- 
state manifold with respect to one particle spin-up along 
the cc-axis at ro, and a second particle spin-up along the 
y-axis at n = (r max ,0i). Figure [TlJ a) shows the net x 
spin distribution and Fig. |ll[ b) shows the net y spin dis- 
tribution in the system at every angular position 9 2 as 
a function of 9\. (The net z correlations are negligible 
and not shown.) As 9\ approaches 90°, two peaks (the 
two remaining particles) emerge in the net x and y spin 
distributions at 9 2 = 180° and 6 2 = 2 70°. As 6 X ap- 
proaches 180°, two different peaks arise at 9 2 = 90° and 
9 2 = 270°. There is a third region, at 9 X = 270° of large 
correlation that mirrors that at 9\ = 90° . Note the inver- 
sion symmetry through the point {9 1 ,9 2 ) = (180°, 180°). 
Figure [12] focuses on the two regions of large correlation, 
when 9\ — 90° and 180°, showing the net x, y, and z spin 
distributions as a function of 9 2 . From these plots we con- 
clude that chiral spin structures exist in the ground-state 
manifold of the interacting four-particle system, but the 



structures cannot be readily characterized by a definite 
winding about any axis. 

We now go on to consider the lowest-energy excitation 
above this ground-state manifold, where we do uncover 
winding textures. 



B. First Excited State 

The first excited state in the four-particle system is 
non-degenerate with quantum numbers L z = 0, S — 0. It 
too has a circularly symmetric single-particle density dis- 
tribution about the origin, with a peak at r max = 48 nm. 

Figure [13] shows the two-point spin correlations 
throughout the plane of the dot with respect to a spin-up 
particle at r = (r max , 0). The distribution of the spin- 
up and spin-down densities are shown. The spin-up den- 
sity shown in Fig. [l~3[a) contains two peaks at 9\ = 90° 
and 9 \ = 270°, both along the ring r max . The proba- 
bility drops to approximately one tenth of its magnitude 
between the peaks, at 9\ = 180°. The spin density goes 
to zero as T\ approaches ro, giving evidence for a Pauli 
vortex at ro. The spin-down density shown in Fig. |13[ b) 
has three peaks, the largest at 9i = 180°, and two smaller 
ones of equal magnitude at 9\ = 90° and 9\ = 270°. In 
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FIG. 13: (Color online.) Two-point spin correlations in the 
first excited state of the four-particle system with respect to 
a spin- up particle at ro = (48 nm, 0) (white cross). Shown 
are the spin-up (a) and spin-down (b) distributions. 
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FIG. 14: (Color online.) (a) Trace of the two-point spin 
distribution along the ring r max = 48 nm with respect to a 
spin-up particle at ro = (48 nm, 0) in the four-particle first 
excited state, (b) Net two-point distribution of S z throughout 
the QD for the same state. 



contrast to the ground state (see Fig. [8]), this manifold 
does not exhibit antiferromagnetic order. 

The two smaller peaks in the spin-down distribution 
are approximately half the magnitude of the large peak 
and approximately equal to the magnitude of the two 
peaks in the spin-up distribution. The net S z distribu- 
tion at 6*i = 90° and 8\ = 270° is therefore approxi- 
mately zero, indicating an in-plane orientation. This is 
further evident in the trace along the ring r max shown 
in Fig. 14 a), and the net spin distribution shown in 
The single peak at 6\ = 180° is composed pri- 



Fig. 14 



marily of a single spin-species. These results are consis- 
tent with two-point spin correlations in Ref. 1101 Of note 
here is that, in the transition from spin-up at &i = 0° to 
spin-down at 8\ — 180°, the spin lies almost completely 
in the plane of the QD. Calculation of the canting an- 
gle along the ring r max as a function of 8\ is plotted in 
Fig. [15] In this case there are two locations where the 
net spin has the greatest canting. Relative to the spin-up 
particle at r = (r max , 0), as 8\ approaches 8o = and 
also as 8\ approaches 180°, the canting angle approaches 
180°. In the regions along r max where the net spin distri- 
bution is almost zero, the canting angle approaches 90°, 
showing that the net spin exists in the x-y plane. How- 
ever, the two-point correlation function is incapable of 
determining the orientation of the spin within the plane. 
We therefore turn to the three-point correlation function 
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FIG. 15: (Color online.) Canting angle a as a function of 6\ 
along the ring r max = 48 nm in the four-particle first excited 
state. The net spin density from the two-point calculation 
along the same ring r max is shown for reference. 
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FIG. 16: (Color online.) Three-point correlations along the 
ring r max = 48 nm in the four-particle first excited state with 
respect to a spin- up particle along the x-axis at ro = (r max , 0) 
and a spin- up particle along the y-axis at ri = (r max , 9i). 
Shown are the spin distributions along r-z = (r max , 62) for 
the net spin along the a;-axis (a) and along the y-axis (b). 



in order to determine the orientation of the in-plane spin 
component as it transitions between spin-up and spin- 
down. 

We calculate three-point spin correlations along the 
ring r max with respect to a spin-up particle along the 
x-axis at ro = (r max , 0), and a spin-up particle along 
the y-axis at ri = (r max ,$i). The net x and net y spin 
distributions are shown in Fig. [16] as a function of 9\ 
and 82- These distributions reveal strong correlations at 
9\ = 90° and 8\ = 270° , consistent with the previous two- 
point correlation. Note the inversion symmetry about 
the point (6>i,0 2 ) = (180°, 180°). The plots further show 
that for B\ = 90°, (270°), the net spin is predominantly 
spin-down along x at 82 = 180°, and predominantly spin- 
down along y at 8 2 = 270° (90°). Focusing attention 
to the case of 8\ = 90°, we plot in Fig. [l7^a) the spin 
distribution along the ring r max given a spin-up particle 
along x at 80 — 0° and a spin-up particle along y at 
81 = 90° . Figure Vf\b) shows the effect of this correlation 
on the net spin distribution throughout the plane of the 
quantum dot. 

These results for the net x and y spin distributions 
are consistent with the z distributions calculated by the 



two-point correlations (see Fig. 14). Those two-point cor- 
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FIG. 17: (Color online.) (a) Trace along the ring r max , reveal- 
ing the net spin distributions in the four-particle first excited 
state given a spin-up particle along x at ro = (r max , 0) and 
a spin-up particle along y at n = (r max , 90°). (b) The cor- 
responding net spin distribution in the plane of the QD. The 
vector field depicts the orientation and the color depicts the 
magnitude. The locations ro and ri are indicated by the blue 
arrows. 



relations revealed a single peak at 180° indicating a spin 
anti-aligned to the one at 9q — along the ring r max . 
The three-point function with respect to a spin-up par- 
ticle along 1 at ro = (r max , 0) and a spin- up particle 
along y at r\ — (r max , 90°) yields a net spin polar- 
ized along the negative x-axis at 9 = 180°, and a net 
spin polarized along the negative y-axis at 6 — 270°, on 
the ring r max . This is evidence of a winding along the 
ring r max . From the symmetry of our correlation oper- 
we can deduce that there are in fact 



ators, Eq. ( 13a 



four orthogonal windings that wind about the origin in 
this manner; two which begin with a particle spin-up 
along the x-axis at ro = (r max , 0) and then differ in the 
direction of their spin polarization along the y-axis at 
i"l = (r max , 90°) (i.e. in the chirality of the rotation), 
and two which begin with a particle spin-down along the 
x-axis at ro = (r max , 0), and again differ in their chirality 
along the ring r max . Thus, we can characterize the first 
excited state of the four-particle droplet as a superposi- 
tion of four different windings, differing by their chirality 
(clockwise and counter clockwise), and by a topological 
charge (±1). (See Ref. SHI for a thorough semiclassical 
description) . 



VI. SUMMARY 

We have computed spin correlations in the two lowest- 
lying states of three-particle and four-particle circular 
two-dimensional QDs to resolve spin textures that ex- 
ist in the system in the presence of strong, long-range 
Coulomb repulsion at zero magnetic field. Our findings 
are summarized in Table [i] From the one-point corre- 
lation, we determine the annular regions of maximum 
spin-density in the QD. As expected, the radial distance 
of this region from the origin depends on the number of 
confined particles, and on the strength of the Coulomb 
repulsion. 

We further compute two-point spin correlation func- 



tions to determine the correlations in the spin field along 
a given direction. Each resulting spin distribution is sym- 
metric through the diameter of the QD on which one par- 
ticle is located. The details of the spin configurations are 
dependent on the unique quantum characteristics of the 
state, particularly the spin and orbital angular momen- 
tum quantum numbers. The two-point correlation calcu- 
lations for each state suggest the presence of a Pauli vor- 
tex at the position of the fixed particle. In addition, these 
results reveal an incipient Wigner phase in the states ex- 
amined above. 

To uncover chiral textures, the three-point spin corre- 
lations are calculated. We plot the spin density with re- 
spect to two particles with mutually perpendicular spins. 
The three-particle states we have investigated exhibit 
splayed textures. The four-particle ground state exhibits 
an antiferromagnetic texture, and the first excited state 
exhibits winding textures. Here, the results indicate that , 
given a particle with spin at a point in the dot, the spin 
field rotates through a plane perpendicular to the origi- 
nal spin orientation with a full 2tt winding as one moves 
along a closed trajectory about the origin. Importantly, 
given the finite size of the system and the full 0(3) spin 
symmetry present, these winding textures are only quasi- 
topological in character. The rotational spin symmetry, 
in particular, allows the spins to continually deform into 
the trivial texture (or, rather, to the ground-state tex- 
ture). In a larger, semiclassical system with at least uni- 
axial anisotropy, the analogous textures will have topo- 
logical character. In fact, in 2DEGs and bulk systems, 
skyr mions and other spin textures have been shown to 

cxi „ t l26l30ti 33l 

In the present small but fully quantum system, incip- 
ient topological structures may manifest themselves if a 
coupling were introduced between the spin field and the 
spatial orientation of the quantum dot itself, analogous 
to the anisotropies typically found in much larger mag- 
netic systems. This could come about through spin-orbit 
coupling. There is evidence that such a coupling may 
break the degeneracy between the spin-winding statesf^l 
Together with the strong Coulomb repulsion present in 
these small quantum systems, the chiral structures that 
emerge should exhibit longer lifetimes and lower decoher- 
ence rates than their more conventional counterparts. 

Improvements in the lifetime of these states could be 
determined, for example, by comparing the relaxation 
rates of the states at and away from the ir deg eneracy 
point by use of single-shot measurements! 51 ^ 52 ! Experi- 
mental methods exist for differentiating between fully 
polarized spin states and correlated spin texture states 
in the ground state of a QD by measuring the excitation 
spectra of the QD as a function of magnetic fieldPStUD 
These methods can be used to determine when a spin 
texture state is present in the system. The issues of spin- 
orbit coupling and state lifetimes will be investigated in 
future work. 
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TABLE I: Summary of textures uncovered in the respective degenerate manifolds. 



Manifold 



Single Particle Densities 



2 Point 



3 Point 



Canting 



3-particle ground state 

3- particle first excited state 

4- particlc ground state 
4-particle excited state 



39 nm 
41 nm 

rmax = 48 nm 
rmax = 48 nm 



triangular lattice 
triangular lattice 
antiferromagnetic order 
in-plane configuration 



splayed 
splayed 
no winding 
in-plane winding 



131°* 
53° 
varies 
varies 

[at peaks) 
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Appendix : Derivation of Product Spin Operators 

In this appendix we discuss the general properties of 
the product-operators used in our above analysis, as 
well as the details associated with the derivation of each 
product-operator. 



1. One-Body Spin Operators 

A general one-body operator is expressed in canonical 
second-quantized form as 



U = J2 Ua P a « 6 /3' 

a /3 



where U a g = ( a \U\(3) is a matrix element of the 
operator p2 and and c a are second-quantized Fermi 
operators, respectively creating and destroying a particle 
in state \a). 

For a system of spin-1/2 fermions, the operator for the 
spin density at position r is given by (h = 1) 



S(r) = ~X)#(r)S"' 



ips'(r), 



(A2) 



where <x = [a x , & v , <r z ) are the Pauli spin matrices and 
ip\{r) and V> s (?") are the field operators. 

Equation (A2) yields the net spin density at point r. 



We arc additionally interested in distinguishing the spin- 
up and spin-down densities along each coordinate-axis. 
We define a general set of spin operators S± a (r) that 
separately determines the spin-up and spin-down densi- 
ties along a = x, y, z at position r. In the S z basis, the 
operator for the ±x spin density, for example, is given by 



S± x (r) 



1 r 

4 



^j(r)^ t (r) + ^j(r)^(r) 
± (#(r)^(r) + #(r$ t (r))] . (A3) 



This operator can be derived from the field opera- 
tors, V-L(r) = $[(r) ± $\(r)]/y/2, and S ±x (r) = 

ip± x (r)ip± x (r)/2. Similarly, the operators along the 
other two orthogonal directions are 



S± y (r) 



and 



S+Ar) 



± i 



(j>l(r)i> t (r) 



7/>j(r)-0 ; (r; 



(A4) 



$\(r)i> t (r), S^ z (r) 



^j(r)V^(r). 



(A5) 

By defining the net spin along an axis to be the difference 
between the spin-up and the spin-down density along 
that same axis (S x (r) = S +X (r ) — S- X (r), etc.), we ob- 
tain the components of Eq. ( |A2[ ) . Upon integrating these 
net spin components over all space, we recover the usual 
expressions S x = (S++S-)/2, S y = (S+-S-)/(2i), and 
(Al) S z = (Sf — Si)/2, where S + and 5_ are S z raising and 



lowering operators, respectively. 

The operators in Eqs. ((A3) [A4| 
which flip spins 



each contain terms 
But since the Hamiltonian, Eq. , con- 



serves spin, and since we include spin quantum numbers 
to classify our states, these terms give zero contribution 
to Eq. (5a) for one-body spin operators. As a conse- 
quence of spin conservation we have for any degenerate 
manifold at zero magnetic field, (S x (r)) — (S y (r)) = 
(S z (r)} = 0. Along the quantization (z) axis, we can 
distinguish between the spin-up and the spin-down den- 
sity at position r, but we cannot distinguish between the 
spin-up and spin-down density along the two orthogonal 
directions: Formally, (S+ X (r)) = (S- X (r)), and similarly 
for the y direction. Consequently, it suffices to calculate 
only the one-point correlations of S± z (r), as in Eq. 



2. Two-Body Spin Operators 

To investigate correlations along a particular axis, the 
two-point correlations functions are required. On phys- 
ical grounds, we require that operators be symmetric in 
their indexes.^ For one-body operators, we require that 
U a p — Uff a . For a two-body operator, we similarly re- 
quire that V a j3 7 s — VpaS-y- In general, a product of 
one-body operators is not an Af-body operator. Let U 
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and V each be a one-body operator. Their product can 
be written as 



At; 



UV = y / (UV) ij c i c j 



ij 



UjtV^clcfac,, (A6) 



ijkl 



where (UV) iS = (i\UV\j), U ik V sl = (i\U \k) (j\V \l). 
The product of two one-body operators is in fact a sum 
of canonical one-body and two-body operators. 

There is no preferred axis along which to calculate our 
two-point correlations because the net spin in each degen- 
erate manifold is zero, but one may exploit the symmetry; 
for example, the correlation between a spin-up particle 
and the remaining spin-up distribution is identical to the 
correlation between a spin-down particle and the remain- 
ing spin-down distribution. In general, the correlations 
between a particle with spin and the remaining particles 
of parallel spin will be the same for any orientation, as 
will the correlations between a particle with spin and the 



remaining particles of antiparallel-spin. Due to spin con- 
servation, two-point correlations between particles with 
perpendicular spin do not provide additional correlation 
information and are therefore not considered. 

The two-point correlations investigated above are 
(S +z (r )S +z ( ri )) and (S +z (r )S- z (n)) in Eq. ([To ) . 
With the condition To 7^ t*i (while still allowing 
ro — > T"i), the one-body term vanishes. As described 
in Sec. |A 1| spin symmetry implies that all four of 
(S +z (r )S± x (r 1 )) and (S +z (r )S± y {r 1 )) are identical 
and equal to ((S +z {r )S +z (n) + S +z (r )S- z (n)))/2. 



3. Three-Body Spin Operators 

To probe for chiral textures, we compute the three- 
point correlation functions. These are a product of three 
one-body operators. In canonical form, the product of 
three one-body operators is the sum of a canonical one- 
body, two-body, and three-body operator, 



uvw = J2( uvw h &j + lj2(( UV ^Wji + U ik (VW) kl + {UW) lk V kl 



ijkl 



+ (UV) 3 iW lk + U jt {VW) ik + {UW) ol V lk )c\cy Cl c k 



jl ^2{ U U V jmW kn + U u V kn W jm + U jm VuW kn 

+ U jm V kn Wii + U kn VuW jm + U kn V jm Wi^ ct44 e nC m C; , (A7) 



ijk 
Iran 



where (UVW)jj = (i\ UVW \ j). Each of the matrix el- symmetry implies that the one-body and two-body pieces 



ements in Eq. (A7| is symmetric under appropriate in- of Eq. (A7| vanish when r , ri, and r 2 are all different 



tcrchange of indexes; for a three-body matrix element 
Viju mn , for example, we require V ijk i mn = V lkj i nm , and 
so on. 

We express our three-point correlation operators in the 



symmetric form shown in Eq. ( A7 1 . As in Sec A 2 spin 



Because our averages are taken with respect to spin- 
conserving states, we need only consider terms in the 
correlation function which themselves conserve spin. We 
can thus write the product operators as they are given in 
Eq. pi. 
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